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ABSTRACT 

Based on three years of deep observations of the Galactic center with the Chandra X-ray Observatory, 
we report the discovery of changes in the intensities and morphologies of two hard X-ray nebulosities. 
The nebulosities are dominated by fluorescent iron emission, and are coincident with molecular clouds. 
The morphological changes are manifest on par sec scales, which requires that these iron features 
are scattered X-rays from a 2 or 3-year-long outburst of a point source (either Sgr A* or an X-ray 
binary) with a luminosity of at least 10 37 erg s _1 . The variability precludes the hypotheses that these 
nebulae either are produced by keV electrons bombarding molecular clouds, or are iron-rich ejecta 
from super novae. Moreover, the morphologies of the reflection nebulae implies that the dense regions 
of the clouds are filamentary, with widths of ~0.3 pc and lengths of « 2 pc. 
Subject headings: Galaxy: center — ISM: clouds — X-rays: ISM 



1. INTRODUCTION 

As a result of the high density of matter and stars in the 
Galactic center, opportunities arise for energetic parti- 
cles, radiation, and interstellar matter to interact in ways 
that occur very rarely in the Galactic disk. One manifes- 
tation of these interactions are iron fluorescence nebulae 
that are visible in the X-ray band, and that are associ- 
ated with the molecular cloud complexes Sgr Bl and B2 
(Murakami et al. 2000, 2001b), Sgr C (Murakami et al. 
2001a), and M 0.11-0.11 (Park et al. 2004). The fluores- 
cence is thought to be produced when low-ionization iron 
in the molecular clouds is bombarded either by ~30 keV 
electrons (Valinia et al. 2000; Yusef-Zadeh, Law, & War- 
die 2002; Bykov 2002) or by -8 keV X-rays (Sunyaev 
& Churazov 1998; Revnivtsev et al. 2004b). However, 
exciting the fluorescent features with electrons would 
require a large, and arguably unreasonable, amount of 
energy in keV electrons (Revnivtsev et al. 2004b), and 
there is currently no X-ray source in the Galactic cen- 
ter bright enough to illuminate these features (Sunyaev 
et al. 1993; Koyamaet al. 1996; Revnivtsev et al. 2004a). 
Therefore, the excitation mechanism for these features is 
uncertain. 

In this paper, we attempt to determine how the iron 
fluorescence features are excited, by searching for varia- 
tions in the morphologies and intensities of the two exam- 
ples with the highest surface brightnesses (Features 1 and 
2 in Park et al. 2004), which are located in the M0. 11- 
0.11 cloud complex. If the fluorescence is excited by par- 
ticles, which would propagate rather slowly (<0.06c for 
30 keV electrons), the features would not vary over the 
course of a few years. Alternatively, if the fluorescence 
is produced by X-rays, the absence of any source that 
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is currently bright enough to produce them implies that 
the X-rays were produced in a transient event. If the 
outburst were short enough, then in a few years the light 
front could propagate across the features, causing their 
morphologies and intensities to change (Sunyaev & Chu- 
razov 1998). Here we report that, in fact, the two fluo- 
rescent features have changed after three years. 

2. OBSERVATIONS 

The features that we studied were identified by Park 
et al. (2004), based on 500 ks of observations of the 20 pc 
around Sgr A* taken with the Advanced CCD Imaging 
Spectrometer imaging array between 2002 May 22 and 
June 3. Since then, 100 ks of data were acquired between 
2004 July 5 and 6, and 200 ks between 2005 July 24 and 
August 1. All of the observations had aim points within 
20" of Sgr A*, although the roll angles for the satellite 
differed in 2002 (76°) and 2004-2005 (280°). Features 1 
and 2 are located near a=266.492°, £=-28.943° (J2000), 
or 6' from Sgr A*, and are covered by all of our obser- 
vations. 5 At this offset, the point-spread function (PSF) 
at 6 keV has a full-width at half-maximum (FWHM) of 
6", which corresponds to a physical scale of 0.2 pc (D=8 
kpc; McNamara et al. 2000). 

We processed the event lists from each observation to 
remove bad pixels, cosmic rays, and background flares 
in the standard manner 6 , using CI AO version 3.3.0 and 
the calibration database version 3.2.2. Next, we identi- 
fied point-like X-ray sources using the wavelet algorithm 
wavdetect. Finally, the absolute astrometric frame was 
registered using X-ray sources with counterparts in the 
2MASS catalog (Skrutskie et al. 2006), making it accu- 
rate to 0"2 . 

To search for temporal evolution, we produced three 
separate images of the flux in the 4-8 keV band, using 
the data from 2002, 2004, and 2005. This limited energy 
range was chosen to reject foreground emission, emission 
from soft plasma, and particle background. In Figure 1, 

5 A third feature was also reported at a:=266.593, 5=-28.920 
in Park et al. (2004), but was not covered by the observations in 
2004-2005, and so we ignore this feature. 

6 http://asc.harvard.edu/ciao/threads/ 
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Fig. 1. — Hard (4-8 keV) X-ray images of the 6'x 6' field around 
the iron fluorescent features (centered at a, S = 266.44156, - 
28.94267 [J2000]), illustrating their evolution as a function of time. 
The solid ellipses denote the locations of the fluorescent iron fea- 
tures from Park et al. (2004). The dashed ellipses in the top panel 
denote the regions used for background extraction. The bar at top 
displays the mapping between flux and color, with red indicating 
the brightest emission, and blue the faintest. 



we display the portion of the field containing Features 
1 and 2, smoothed by a Gaussian with a FWHM of 6". 
Having smoothed the images on this scale, changes in 
the PSF shape are not significant, as can be seen by 
examining point sources in the image. Over three years, 
the morphologies of Features 1 and 2 change noticeably. 
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Fig. 2. — Significance maps of the difference in flux received 
in the 4-8 keV band. Red indicates an increase in intensity. At 
the locations of the fluorescent features, the variations in flux are 
spatially coherent, indicating that the diffuse features have changed 
in shape and intensity. The changes are more significant in the 
2005-2002 difference map, because the exposure was twice that 
in 2004. The bar at the top of the figure provides the mapping 
between significance and color. 



These variations can be seen more clearly in difference 
images, which we have normalized to the uncertainty in 
the convolved flux maps and displayed in Figure 2. To 
generate these, we first calculated the one-sided, upper 
la Poisson uncertainty on the counts in each pixel of each 
image, following Gehrels (1986). We divided this by the 
corresponding exposure maps to obtain the uncertainties 
on the flux, <jf^. The uncertainties in each convolved 
image ((Jc,i) were then the square root of the convolution 
of the square of the uncertainties with the square of the 

Gaussian smoothing function (#), ac,i = (cr'j * g 2 )]^ 2 • 
The uncertainties in each pair of images were then added 
in quadrature to obtain those of the difference images. 
The difference images were divided by the uncertainty 
maps to obtain the significance maps. 

In Figure 2, red features were defined to be brighter 
in the later images; a clear example is the transient 
source XMM J174554.4-285456 (Porquet et al. 2005; 
Muno et al. 2005), which appeared in 2004 just north 
of the center of the image. We find significant variations 
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TABLE 1 

Spectra of the Iron Fluorescence Features 
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2002 2004-2005 



Feature 2 
2002 2004-2005 
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Energy (keV) 
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Energy (keV) 



Fig. 3. — Spectra of fluorescent Features 1 and 2. The top panels 
display the spectra in units of detector counts (crosses). The solid 
line is the best-fit model of an absorbed power law and two narrow 
lines at 6.4 and 7.0 keV. The bottom panels display the residuals, 
in units of the statistical uncertainty in the counts in each bin. The 
dominant feature in each spectrum is the iron line at 6.4 keV. 



in several places in the field. Some of the most promi- 
nent variations are at the locations of the two fluorescent 
features with the highest surface brightness. The main 
ridge that composed Feature 1 in 2002 faded at the So- 
level by 2005, while the northern part of the feature on 
either side of the ridge brightened by 4<j. For Feature 2, 
we find that the southern portion faded by 3<r between 
2002 and 2005, whereas the eastern portion shifted far- 
ther east and became more extended at the 4-5<r level. 
We also find 3-4a increases in the flux to the east of Fea- 
ture 2 and the south of Feature 1. The pervasive nature 
of these variations is not surprising, becuase then entire 
field contains molecular gas (Tsuboi et al. 1999) and ex- 
hibits diffuse iron Ka emission (Muno et al. 2004). 

Next, we extracted spectra of the diffuse features. 
The source spectra were extracted from the solid el- 
lipses in Figure 1 (top panel) and the background from 
the dotted ellipses above and below the features. We 
excluded events that fell within circles circumscribing 
the 92% encircled-energy contours of the PSFs for each 
point source. The background-subtracted spectra were 
grouped so that each bin had a signal-to-noise ratio of 
>5. Spectra for all of the data are displayed in Figure 3. 

We modeled the spectra from the 2002 and 2004-2005 
epochs jointly using XSPEC, as a power-law continuum 
with photon index T and normalization 7V r , two nar- 
row iron lines at 6.40 (Fe Ka) and 7.06 keV (Fe K/3). 
We assumed that this emission was absorbed by the in- 
terstellar medium (-/Vhjsm)? plus a component of the 
absorption intrinsic to the associated molecular clouds 
(^h, cloud) that covered a variable fraction of the emit- 
ting region (/cloud)- The use of such a "partial covering" 
factor (/cloud) was motivated by the assumption that the 
emission from the clouds would take place at a range 
of optical depths. We also assumed that T and 7Vh,ism 
were identical for both features during both epochs. In 
Table 1, we list the best-fit parameter values, iron equiv- 
alent widths, and fluxes. 

As found by Park et al. (2004), both features have ex- 
tremely strong Fe Ka emission, with equivalent widths of 
wl keV. The Fe Kf3 emission is wlO-15% of the strength 
of the Ka line, as is expected for fluorescent emission 
(e.g., Murakami et al. 2000). Both features are heavily 
absorbed, so most of the observed flux is received be- 
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Note. — The reduced chi-squared for the joint fit was 211/176. 
The units are as follows: A^h,ism and A^h, cloud are 10 22 cm -2 ; Nr 
is 10 -4 ph cm -2 s _1 keV -1 at 1 keV; F Ka and F K/3 are 10 -6 ph 
cm -2 s _1 ; EW K a and are eV; F x is the observed 4-8 keV 

flux in units of 10 -13 erg cm -2 s _1 ; and F^ is the absorption- 
corrected 4-8 keV flux in units of 10 -13 erg cm -2 s _1 . Uncertain- 
ties are lcr, for 1 degree of freedom for the parameters, and for 7 
degrees of freedom for the fluxes and equivalent widths. 

a The interstellar absorption and photon index of the continuum 
were assumed to be identical for each spectrum. 

b For Feature 1, / c i ud was poorly-constrained, so we fixed the 
parameter to its best-fit value of 1. 



tween 4 and 8 keV. We find that for Feature 1, despite 
the noticeable changes in the morphology, the 4-8 keV 
flux and the iron flux in the extraction region remained 
constant between 2002 and 2005. For Feature 2, the in- 
ferred iron flux remained constant, but the observed 4-8 
keV flux declined by 26±7%. 

3. DISCUSSION 

We have found that two fluorescent iron features lo- 
cated 6' (14 pc in projection) from Sgr A* varied in mor- 
phology and intensity between 2002 and 2005 (Figs. 1 
and 2, and Tab. 1). These changes in the fluorescence 
features require that the illuminating source varied on 
a time scale shorter than 3 years. Moreover, the fea- 
tures are ^0.3' wide and w V long, which corresponds 
to projected physical sizes of ^0.7 x 2 pc. The smaller 
of these dimensions constrains the minimum duration of 
the outburst to have been two years. The fact that two 
features with a projected spearation of ~12 light-years 
varied over the same time interval implies that the illu- 
minating source and the target clouds have separations 
with a component along our line of sight. 

Such rapid changes over the course of a few years would 
not be produced by steady, slow-moving sources such 
as ejecta from supernovae (Bykov 2002). These shifting 
iron fluorescence features also would not be produced by 
~30 keV electrons, which would have velocities of only 
0.06c (and slower diffusion speeds), and which would pre- 
sumably be accelerated in shocks that persist for much 
longer than a few years. The most reasonable explana- 
tion for these variable fluorescent features is that they 
are scattered emission from the outburst of a transient 
hard X-ray source, such as an X-ray binary or the super- 
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massive black hole Sgr A*. 

The luminosity of the source illuminating these reflec- 
tion nebulae can be estimated from the optical depths 
of the scattering media coupled with the fluxes in the 
iron lines. The luminosity in the fluorescent Fe Ka line 
is about 30% of the flux absorbed by the photoelectric 
edge, and so the observed flux can be computed from 

Fk - 1'^ All- e -r F e[E/Eo}-*] dE (1) 

tKa ~ AttD* J Eo 4^ 2 I 6 J d ' [ ' 

where D is the distance from the observer to the fluo- 
rescence region, Cx{E) is the specific luminosity of the 
illuminating source as a function of energy (units of pho- 
tons s _1 keV -1 ), d is the distance between the source 
and the fluorescing cloud, A is the area of the cloud pro- 
jected toward the source, and rp e is the optical depth of 
the iron photoelectric edge at an energy Eq = 7.1 keV. 

This equation can be simplified given the parameters 
of our nebulae. The optical depth is given by t p e = 
A^H^Fe^Fe, where <TF e = 3.5 x 10 ~ 20 cm 2 is the cross- 
section for photoelectric absorption at 7.1 keV, and 5p e = 
3 x 10 -5 is the fractional abundance of iron for solar 
abundances. If we take iV H < 6 x 10 23 cm -2 (Tab. 1), 
the optical depth to iron absorption is rp e < 0.8, and 
1 — e~ T « r with an error of <30%. Next, we take D=S 
kpc, define dio to be normalized to 10 pc, re-write A as 
a solid angle normalized to 1 arcmin 2 at 8 kpc, and 
write t Fe in terms of the column depth A^h,23 in units 
of 10 23 cm -2 and the abundance of iron relative to the 
solar value ap e - Finally, we assume that the spectrum of 
the illuminating source is a T=1.8 power law, Cx(E) oc 
E~ 1,8 . We then can integrate eq. 1 over energy, and write 
the result in terms of the 2-8 keV luminosity, normalized 
to 10 38 erg s _1 , £x,38- With these substitutions, we find: 

F Ka = 5 x 10~ 5 L x ,38 Nu,23 «Fe ^i ph cm -2 s _1 . 

(2) 

If Sgr A* is the source of the photons that produced the 
reflection features, and that the line-of-sight separation 
between it and the cloud is equal to the projected separa- 
tion (14 pc), d « 20 pc and an outburst with Lx ~ 10 38 
erg s _1 would have been required. 7 Such an outburst 
would have occurred 60 years ago, before the advent of 
X-ray astronomy, and therefore would not have been ob- 
served. The lack of similar reflection nebula at negative 
Galactic longitudes (Park et al. 2004) implies that those 
clouds do not lie within the light front of the outburst. 

Alternatively, if the bright transient that was in out- 
burst in 2004 produced the fluorescent features, d « 7 
pc and Lx ~ 10 3T erg s _1 . Although this is 100 times 
brighter than the outbursts in 2002 and 2004 (Porquet 
et al. 2005; Muno et al. 2005), this outburst would have 
occurred 20 years ago, before the advent of both Chandra 
and the wide-field X-ray monitors on BeppoSAX and the 
Rossi X-ray Timing Explorer, so it too could have been 

7 This luminosity is larger than was reported in Park et al. 
(2004), because a mathematical error was made in computing the 
value in the earlier paper. 

8 See also http://xte.mit.edu/ASM_lc.html. 



missed. However, it is uncommon for the outburst of an 
X-ray binary to persist at this luminosity for 2-3 years 
(Chen et al. 1997) 8 , so we prefer the hypothesis that the 
source of the outburst was Sgr A*. 

The time evolution of a reflection nebula has been ex- 
plored in detail by Sunyaev & Churazov (1998), assum- 
ing target material with a spherical geometry. However, 
the observed changes in the morphology in Figure 1 are 
significantly more complex than they considered. For in- 
stance, in Feature 1 the vertical filament ~2 parsec long 
fades over the course of 3 years. The fact that the light- 
travel time along the feature is longer than the interval 
over which it fades implies that the light front propa- 
gated perpendicular to the filament, and had width that 
is similar to that of the filament (0.3 pc). The fact that 
the whole filament dimmed, with very little brightening 
adjacent to most of its length (aside from knots in the 
north on either side of the filament) implies that the X- 
rays were no longer propagating through dense material 
in 2004 and 2005. Therefore, the morphology of the iron 
features traces the structure of the clouds, revealing that 
they contain filamentary features a few tenths of a parsec 
wide and a couple parsecs long. 

4. CONCLUSIONS 

The discovery of changes in the intensity and morphol- 
ogy of fluorescent iron nebulae near the Galactic center 
demonstrates that they are reflection nebulae. The most 
likely source of the reflected X-rays is an outburst of 
Sgr A* with a luminosity of at least 10 38 erg s _1 , which 
would have occurred about 50 years ago and lasted 2-3 
years. If a similar outburst had occurred since then, it 
would have had a good chance of being detected by an 
X-ray instrument, or would have produced similar fluo- 
rescent features in other, nearby clouds. Since no other, 
long outburst of Sgr A* has been reported, and no flu- 
orescent features have been identified closer to Sgr A* 
(Park et al. 2004), the duty cycle of bright, year- long 
outbursts from Sgr A* must be <5%. 

These fluorescent features also reveal filamentary 
structures within the dense regions of the molecular 
clouds. Inter ferometric millimeter observations would 
resolve clouds on a similar spatial scale, and addition- 
ally provide velocity information. Combining the cloud 
velocities and the three-dimensional spatial information 
afforded by the evolution of the reflection nebulae would 
allow one to construct a physical model of these turbulent 
molecular clouds. This is a necessary step toward mod- 
eling star formation in the hot, turbulent environment of 
the Galactic center. 
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